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A  Comparison  of  Three  Methods  for  Computing  the 
Added  Mass  of  Ship  Sections 


Abstract 

Three  methods  that  have  been  developed  for  computing  added-mass 
coefficients  of  two-dimensional  forms,  oscillating  horizontally  or  verti¬ 
cally  at  a  free  surface  are  compared  by  application  to  a  set  of  four  cargo- 
ship  sections.  Computer  programs  are  listed  for  two  of  these  methods.  The 
method  employing  conformal  mapping  is  recommended  as  the  best  of  the  three. 

Introduction 

In  this  work,  a  comparison  of  three  methods  previously  developed 
by  Landweber  and  Macagno,  [l,  2,  3]  for  computing  added-mass  coefficients 
of  two-dimensional  forms  oscillating  at  a  free  surface,  will  be  presented. 

In  the  first  of  these  papers ,  general  expressions  for  added  mass 
for  horizontal  and  vertical  oscillations  were  derived,  but  were  applied 
only  to  a  two-parameter  family,  the  so-called  Lewis  forms.  In  order  to 
extend  the  family  of  Lewis  forms,  a  three-parameter  family  was  treated  in 
reference  [2].  In  practice  the  two-parameter  method  has  been  found  more 
convenient  and  has  been  more  generally  used.  Since  both  of  these  methods 
are  indirect  and  approximate,  a  procedure  employing  conformal  mapping, 
which  yields  the  actual  added-mass  coefficients  of  a  given  section,  was 
developed  and  reported  in  reference  [3]. 

In  order  to  encourage  the  wider  use  of  the  three-parameter  method 
which  is  intrinsically  more  accurate  than  the  two-parameter,  and  to  clarify 
the  application  of  the  conformal  mapping  procedure,  it  seemed  necessary  to 
prepare  complete  computing  programs  for  the  latter  two  methods. 

The  purpose  of  this  work  is  then  two-fold;  first  to  present  these 
computer  programs  and  second,  by  applying  them  to  a  family  of  ship  sections 
to  determine  how  well  the  results  by  the  three  methods  agree,  and  which  to 
recommend  for  future  use. 


Methods  of  Determining  Added  Masses 
(a)  The  two-parameter  technique 


In  the  first  paper  [l],  the  general  theory  was  applied  to  obtain 
expressions  for  the  added-mass  coefficients  of  arbitrary  two-dimensional 
forms,  on  the  assumption  that  the  conformal  transformation  of  the  exterior 
of  the  closed  form  into  the  exterior  of  the  unit  circle  in  the  ^-plane  is 
known  and  given  by  the  expression 


where 


real 


(1) 


z=x+iy  ,  £=£+in 

A  particular  family  of  ship-like  sections,  the  Lewis  forms,  were  then 
derived  from  the  unit  circle  by  this  transformation  for  the  case  where  a^ 
and  a^  are  the  only  nonzero  coefficients.  In  parametric  form,  the  equa¬ 
tions  of  these  forms  are 


x  =  (l  +  a  )  cos  0  +  a0  cos  30  ) 

13  (2) 

y  =  (1  -  a1)  cos  0  -  sin  30  J 

If  b  denotes  the  half-beam  of  the  section  at  the  water  line,  and  H  the 
draft  at  the  keel  line,  the  coefficients  a^  and  a^  are  related  to  b 
and  H  by  the  expressions 


b  =  1  +  a^  +  a^ 


H  =  1  -  a^  +  a^ 


or ,  solving  for 


and  a 


3’ 


,  1-X  .1+1  /  \ 

a.±  -  b  -Tj-  a3  =  b  -j-  (3J 

where  X  =  —  • 
b 

If  H  and  b  are  given,  the  sectional  area  of  the  form  is  ob- 


S  =  |[-  2  +  3b (l  +  X)  -  b2(l  +  A  +  X2)] 


tained  from 


-3- 


and  the  section-area  coefficient,  defined  by 


is  given  in  terms  of  a  parameter  a  =  2/b,  by 

a  =  qy  [-  a2  +  3a(l  +  X)  -  2(l  +  X  +  X2)] 


If  the  shape  characteristics  X  and  o  are  given,  the  corresponding 
values  of  b,  H,  a^,  a  and  a  can  be  obtained,  but  not  all  values  of  a 
give  possible  useful  forms;  i.e.,  a  must  lie  within  a  certain  range.  Once 
a  is  found,  the  coefficients  for  vertical  vibrations,  C  ,  and  for  hori- 

7  Y 

zontal  vibrations  C„  are  given  by 

n 

C  =  1  +  (1  +  X  -  a) ( X  -  a) 

°H  =  ^[1  +  3X^  (1  +  X  "  “)2] 


A  series  of  Lewis  forms  and  corresponding  curves  for  the  added-mass  coef¬ 
ficients  are  given  in  this  paper.  Therefore,  given  X  and  a  for  any  ship- 

section  that  resembles  a  Lewis  form,  the  values  of  C  and  CTT  can  be 

v  H 

computed  from  the  formulas  (L)  to  an  accuracy  that  depends  on  how  closely 
the  Lewis  form  represents  the  given  ship  section. 


As  was  shown  by  Prohaska  [L],  the  two  parameters  draft-beam  ratio 
X  and  section-area  coefficient  a  are  insufficient  to  define  an  added  mass 


and  since  Lewis  forms  cannot  represent  ship  sections  with  area  coefficients 
close  to  unity,  a  new  method  was  developed  [2].  In  this  paper,  a  more 
general  three-parameter  family  of  forms  was  derived,  for  which  the  added- 
mass  coefficients  were  determined  and  presented  as  a  series  of  curves. 

This  new  family  of  forms  is  a  natural  extension  of  the  Lewis  forms;  to  the 
two  parameters  that  were  previously  used,  a  third  parameter  n  was  added, 
the  ratio  of  the  radius  of  gyration  about  the  transverse  axis  in  the  free 


surface  to  the  draft. 
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In  this  case,  the  family  of  forms  was  defined  from  the  unit  cir¬ 
cle  in  the  C-plane  by  the  transformation 


z 


f5 

C5 


a^,  a^,  a,_  real 


The  equations  of  the  form  are  now 

x  =  (1  +  a1)  cos  0  +  a^  cos  30  +  cos  50 
y  =  (l  -  a1)  sin  0  -  a^  sin  30  -  a^  sin  50 


(5) 


If,  as  before,  b  denotes  the  half  beam  of  the  form  and  H  the  draft  at 
the  keel  line,  from  equations  (5)  we  obtain 


b  =  1  +  a1  +  a^  +  a^ 

H  =  1  -  a  +  a  +  a,- 
13  5 

The  parameters ,  in  terms  of  b ,  H  and  S ,  are  given  by 


(6) 


2  H  _S_ 

"  h  ’  A  "  b  *  a  -  2bH 


(T) 


where  I  is  the  moment  of  inertia  about  the  x-axis,  given  by 

H  2  * 

1  =  2/  xy2dy  =  -  rr  /  y3  dx 

0  0 


(8) 


in  which  the  second  form  is  obtained  from  the  first  by  integration  by  parts . 

A  quadratic  expression  relates  the  parameters  a,  X,  a  and  a,.. 


a  =  {-  a2(l  +  3a,-2 )  +  a[3(l  +  x)  +  a^(l  -  x)]  - 

2(1  +  X  +  X2)}  (9) 


Solving  for  a^  and 


from  equations  (6),  we  obtain 
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and  from  (7),  (8)  and  (10)  we  obtain  for  n 

k 

n  =  j_28 +  ^2&3  +  2^^a32  ~  a3a5  +  a52^  ~  a32a5^  ~ 

-  3^  -  233a3  -  232(4a32  -  5a3a^  +  6a52)  + 

+  123a32a^  -  (3a32  +  a,_2)(a32  +  5a,-2)} 


(10) 


(11) 


where  3  =  1  -  a  . 

When  the  parameters  X,  a  and  a^  are  given,  the  corresponding 
values  of  b,  H,  a^,  a3  and  n  can  be  computed  from  (9),  (10)  and  (6). 
Finally,  expressions  for  and  C^,  in  terms  of  those  parameters  are 

given  by 


C  =  77  {3<*2a  2  -  a(a  -  X  +  l)a  +  2[(X  -  a)2  +  X  -  a  +  l]  } 
v  2  5  5 


,  .  aaq 

CH  =  ^  {[1  +  ^(1  +  X  -  a)2]  +  [5  +  bOX  +  a(l+0  +  156a5)]} 


(12) 


To  use  this  technique  for  computing  and  C^,  the  following 

procedure  may  be  applied: 

l)  Compute  I  from  (8).  If  the  ship  form  is  given  in  polar  co¬ 
ordinates  (r,  4>),  where  x  =  r(<|>  )cos<}> ,  y  =  r(<j>)sin4>  then,  since  dy  = 
cLx* 

=  (r  cos<|>  +  -rr  sin<j>)d<f>,  (8)  yields,  after  simplification, 
d<p 

n  h 

I  =  —  /  r  sin24>  d<j) 

2  0 


Then  obtain  n  from 


n 
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2)  Assume  a  set  of  values  of  in  the  range  permitted  by  the 
condition  given  by  [2]  and  compute  corresponding  values  of  a  from  (9)  and 
of  a^  and  a^  from  (10). 

3)  Obtain  values  of  n  from  (ll)  for  each  set  of  values  of  a^ ,  a^, 
a^,  a  and  X. 

4)  By  interpolation  among  the  values  of  n  from  step  3,  obtain  the 
value  of  a^  corresponding  to  the  actual  value  of  n  for  the  ship  form, 
obtained  in  step  1. 

5)  Apply  the  a^  computed  in  step  4  to  obtain  the  corresponding 

values  of  a,  an ,  a  from  (9)  and  (10),  and  finally  C  and  C  from 
13  V  n 

(12). 

A  Fortran  program  following  these  steps  is  listed  in  the  appendix. 
In  this  program,  the  input  data  are  the  weighing  factors  for  Simpson's  rule, 
A,  the  polar  radii  r,  the  draft-beam  ratio,  X,  the  section-area  coeffi¬ 
cient  a  and  the  assumed  value  of  a^.  An  IBM7044  takes  forty  seconds  for 
compilation  and  loading  this  program  and  less  than  a  second  for  computing 

C  and  CTT  for  each  form, 
v  H 

(c )  The  conformal -mapping  technique 

In  the  first  two  methods,  it  is  assumed  that  a  ship-section  hav¬ 
ing  the  same  principal  geometric  characteristics  as  a  member  of  one  of  the 
particular  families  (the  Lewis  forms  or*  the  more  general  three-parameters 
forms)  will  have  the  same  added-mass  coefficients  as  that  member.  But,  as 
this  assumption  is  not  exact,  a  third  method  was  derived  [3].  In  this  me¬ 
thod,  the  coefficients  of  the  mapping  function  (l)  are  obtained  directly 
by  conformal  mapping  of  the  given  section  in  the  z-plane  into  a  unit  circle 
in  the  4-plane.  With  the  availability  of  high-speed  computers  it  is  now 
possible  to  perform  conformal  transformations  for  functions  with  a  rather 
large  number  of  coefficients. 

Following  a  method  proposed  by  Bieberbach  for  conformal  mapping 
of  the  interior  of  regions  of  arbitrary  form,  corresponding  results  were 
derived  by  Landweber  and  Macagno  [3]  for  exterior  regions.  To  improve  the 
accuracy,  the  conformal  transformation  is  applied  twice,  because  it  was 
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found  that  a  single  direct  application  of  the  Bieherhach  method  gave  re¬ 
sults  for  the  added  masses  that  were  not  sufficiently  accurate. 


The  expressions  obtained  for  these  added-mass  coefficients  are 


as  follows : 


C  =  A  [2(1  +  a  )  -  -  ] 

V  1  7T 


(13) 


C„  = 


1 6 


H  tt2H2  L  11V  1 


[  0L,  n  (a.,  -  l)  2  +  2a12a^  +  2aq3(a1  -  l)aq  +  •  •  ■  • 


+  C22a32  +  2a23a3a5  ■"••••  1 


or,  with  a  '  =  a^  -  1,  a m'  =  am>  m  -  2,  3,  4,  .. 


l6 


r  r  1 

Crr  =  — rtrpr  )  }  a  a  a 
H  it2H2  l  l  rs  r 


»  r 

S 


(lU) 


where  a  is  defined  by  the  finite  series 
rs 

(2r  -  l)(2s  -  l)  r  1  .  1 .  .  1  + 

ars  "  “  2(s  -  r)(2s  +  2r  -  1)  L  2(2r  -  IT  2r  +  1  2r  +  3 

...  +  2s  _  3  +  2(2s  -  l)  r  <  S 

As  it  is  seen  from  (13)  and  (lU),  only  one  coefficient,  a1 ,  is  required 
for  calculating  C^,  but  all  the  coefficients  of  the  mapping  function  are 
required  for  the  calculation  of  C^. 

A  Fortran  program  written  for  this  method  requires  as  input  data 

the  weighing  factors  for  Simpson's  rule.  A,  and  the  polar  radii  of  a  ship 

section.  In  the  7044  computer  it  takes  fifty-four  seconds  for  compilation 

and  loading  the  program,  and  ten  seconds  of  execution  time  to  give  the 

coefficients  C  and  CTT  for  each  form, 
v  H 

Application  to  Selected  Ship  Sections 

In  order  to  compare  the  three  methods,  four  transverse  sections 
of  a  cargo  ship  were  used,  applying  data  furnished  by  the  David  Taylor 
Model  Basin.  Their  characteristics  are  given  in  Table  1,  and  the  polar 
radii  in  Table  2. 
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TABLE  1 

Geometric  Properties  and  Characteristics 
of  Four  Sections  of  a  Cargo  Ship 


Frame  No. 

Half -Be am 

Draft 

S 

A 

0 

n 

6l 

31.12 

19.27 

938.52 

.6192 

.7825 

.40896 

91 

39.00 

19-99 

1430.22 

.5126 

.9171 

•55586 

170 

35.64 

21.89 

1111 . 4l 

.6143 

.7121 

.34245 

203 

11.30 

22.69 

290.09 

2.0080 

.5657 

.24517 

Fortran  programs  written  for  the  IBM  7044  for  the  three- 
parameters  and  conformal -mapping  methods  were  applied.  The  values  obtained 
for  S,  A,  a  and  n  are  also  listed  in  Table  1.  Values  of  the  coefficients 
of  the  transformations  for  the  first  two  methods  a^,  a^  and  a^ ,  a^ ,  a^ 
respectively  are  listed  in  Table  3. 

Introducing  these  coefficients  in  the  equations  of  the  forms , 
given  by  (2)  and  (5)  respectively,  yields  the  approximating  ship  sections. 
The  given  ship  sections  as  well  as  the  approximations  are  graphed  in  Figs, 
la ,  b ,  c ,  d . 

In  Table  4  the  values  of  C^  and  C^  computed  by  the  three 
methods  are  listed. 

Discussion  and  Conclusions 

The  conformal -mapping  method  has  the  advantage  of  making  use  of 
finding  directly  the  coefficients  of  the  conformal  transformation  of  an 
arbitrary  ship  form.  In  the  other  methods  it  is  assumed  that  the  ship 
form  is  obtained  by  a  transformation  with  two  or  three  fixed  coefficients. 

Evaluation  of  the  results  presented  in  Table  4  on  the  basis  that 
the  conformal-mapping  method  gives  the  most  nearly  exact  values  (as  has 
been  shown  for  sections  where  exact  values  are  known,  see  Table  3,  ref. 
[3]),  shows  that  for  C  ,  the  three-parameter  is  preferable  to  the  two- 
parameter  method.  For  C^,  however,  one  method  is  not  consistently 
superior  to  the  other,  and  the  deviations  of  either  from  the  conformal 
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mapping  values  may  be  quite  large. 

Consequently,  it  is  recommended  that  the  conformal -mapping 
procedure  he  used  for  obtaining  the  added-mass  coefficients  of  ship  sec¬ 
tions  . 
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TABLE  2 


Radii  of  the  Ship-Sections  for  0  <  <f>  <  90° 
Frame  6l 

Radius  Radius 

Deft. _  ft  _  Deft. _ ft 


0 

31.12 

1 

30.92 

2 

30.76 

3 

30.60 

4 

30.38 

5 

30.24 

6 

30.06 

7 

29.90 

8 

29.78 

9 

29.62 

10 

29.48 

11 

29.30 

12 

29.14 

13 

29.02 

l4 

28.88 

15 

28.70 

16 

28.54 

IT 

28.38 

18 

28.24 

19 

28.10 

20 

27.94 

21 

27.76 

22 

27.58 

23 

27.42 

24 

27.28 

25 

27.10 

26 

26.94 

27 

26.78 

28 

26.60 

29 

26.42 

30 

26.24 

31 

26.08 

32 

25.90 

33 

25.70 

34 

25.54 

35 

25.36 

36 

25.16 

37 

25.02 

38 

24.84 

39 

24.68 

40 

24.52 

4l 

24.34 

42 

24.16 

43 

24.00 

44 

23.86 

45 

23.68 

46 

23.50 

47 

23.38 

48 

23.18 

49 

23.04 

50 

22.88 

51 

22.72 

52 

22.58 

53 

22.42 

54 

22.26 

55 

22.12 

56 

22.00 

57 

21.84 

58 

21.70 

59 

21.56 

60 

21.44 

6l 

21.30 

62 

21.20 

63 

21.08 

64 

20.96 

65 

20.82 

66 

20.74 

67 

20.62 

68 

20.54 

69 

20.42 

70 

20.32 

71 

20.22 

72 

20. 14 

73 

20.06 

74 

20.00 

75 

19.92 

76 

19-84 

77 

19.78 

78 

19.72 

79 

19.66 

80 

19.62 

81 

19.54 

82 

19.50 

83 

19.44 

84 

19.40 

85 

19.36 

86 

19.34 

87 

19.32 

88 

19.30 

89 

19.28 

90 

19.27 

TABLE  2  Continued 


Frame  91 


Deg. 

Radius 

ft 

Deg. 

Radius 

ft 

0 

39.00 

46 

27.44 

1 

38.98 

47 

27.16 

2 

38.98 

48 

26.76 

3  • 

38.96 

49 

26.40 

it 

38.96 

50 

25.98 

5 

38. 9^ 

51 

25.64 

6 

38.94 

52 

25.24 

7 

38.92 

53 

24.96 

8 

38.92 

54 

24.62 

9 

38.90 

55 

24.32 

10 

38.90 

56 

24.02 

11 

38.86 

57 

23.76 

12 

38.82 

58 

23.50 

13 

38.74 

59 

23.24 

lit 

38.68 

60 

23.04 

15 

38.60 

6l 

22.80 

16 

38.54 

62 

22.60 

17 

38.42 

63 

22.4  0 

18 

38.32 

64 

22.22 

19 

38.20 

65 

22.04 

20 

38.04 

66 

21.88 

21 

37.90 

67 

21.72 

22 

37-72 

68 

21.56 

23 

37.54 

69 

21.40 

24 

37.32 

70 

21.28 

25 

37-06 

71 

21.14 

26 

36.82 

72 

21.02 

27 

36.52 

73 

20.90 

28 

36.22 

74 

20.82 

29 

35.90 

75 

20.72 

30 

35-56 

76 

20.62 

31 

35.12 

77 

20.54 

32 

34.74 

78 

20.46 

33 

34.28 

79 

20.40 

34 

33.78 

80 

20.34 

35 

33.36 

81 

20.30 

36 

32.78 

82 

20.22 

37 

32.36 

83 

20.18 

38 

31.78 

84 

20.14 

39 

31.24 

85 

20.10 

40 

30.76 

86 

20.08 

4l 

30.20 

87 

20. 04 

42 

29.62 

88 

20.02 

43 

29.10 

89 

20.00 

44 

28.56 

90 

19.99 

45 

28.10 
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TABLE  2  Continued 
Frame  170 


Deg. 

Radius 

ft 

Deg. 

Radius 

ft 

■  i  ■■  i  -»  ■«'  ■  ■ 

0 

35.64 

46 

25.04 

1 

35.32 

47 

24.90 

2 

34.98 

48 

24.80 

3 

34.64 

49 

24.64 

4 

34.32 

50 

24.56 

5 

34.00 

51 

24.40 

6 

33.68 

52 

24.30 

7 

33.36 

53 

24.16 

8 

33.04 

54 

24.04 

9 

32.72 

55 

23.90 

10 

32.40 

56 

23.80 

11 

32.14 

57 

23.70 

12 

31.86 

58 

23.56 

13 

31.60 

59 

23.48 

14 

31.32 

60 

23.38 

15 

31.06 

6l 

23.30 

16 

30.78 

62 

23.22 

17 

30.50 

63 

23.14 

18 

30.28 

64 

23.04 

19 

30.02 

65 

22.96 

20 

29.76 

66 

22.90 

21 

29.54 

67 

22.80 

22 

29.28 

68 

22.74 

23 

29.08 

69 

22.68 

24 

28.86 

70 

22.62 

25 

28.64 

71 

22.54 

26 

28.42 

72 

22.50 

27 

28.24 

73 

22.44 

28 

28.02 

74 

22.40 

29 

27.82 

75 

22.34 

30 

27.64 

76 

22.28 

31 

27.42 

77 

22.24 

32 

27.26 

78 

22.18 

33 

27.10 

79 

22.16 

34 

26.88 

80 

22.12 

35 

26.70 

81 

22.08 

36 

26.52 

82 

22.08 

37 

26.40 

83 

22.04 

38 

26.20 

84 

22.00 

39 

26.08 

85 

21.96 

40 

25.96 

86 

21.94 

4l 

25.76 

87 

21.93 

42 

25.62 

88 

21.92 

43 

25.46 

89 

21.91 

44 

25.32 

90 

21.90 

45 

25.18 
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TABLE  2  Continued 


Frame  203 


Deg. 

Radius 

ft 

Deg. 

Radius 

ft 

0 

11.30 

46 

10.22 

1 

11.18 

47 

10.34 

2 

11.00 

48 

10.44 

3 

10.86 

49 

10.56 

4 

10.72 

50 

10.68 

5 

10.58 

51 

10.82 

6 

10.1+6 

52 

10.98 

7 

10.31+ 

53 

11.10 

8 

10.24 

54 

11.28 

9 

10.12 

55 

11.42 

10 

10.00 

56 

11.60 

11 

9.94 

57 

11.80 

12 

9.86 

58 

12.04 

13 

9.80 

59 

12.28 

l4 

9.72 

60 

12.50 

15 

9.68 

6l 

12.76 

16 

9.62 

62 

13.04 

17 

9.58 

63 

13.36 

18 

9*54 

64 

13.68 

19 

9.50 

65 

14.00 

20 

9.48 

66 

14.34 

21 

9.42 

67 

14.72 

22 

9.40 

68 

15.16 

23 

9.38 

69 

15.64 

2k 

9.34 

70 

16.08 

25 

9.32 

71 

16.54 

26 

9.32 

72 

17.02 

27 

9.32 

73 

17.52 

28 

9.32 

74 

18.00 

29 

9-32 

75 

18.50 

30 

9.32 

76 

19.04 

31 

9.32 

77 

19.60 

32 

9.36 

78 

20.08 

33 

9.40 

79 

20.58 

34 

9.42 

80 

20.96 

35 

9.44 

81 

21.36 

36 

9.52 

82 

21.74 

37 

9.54 

83 

22.06 

38 

9.60 

84 

22.36 

39 

9.66 

85 

22.60 

ko 

9.72 

86 

22.76 

Ul 

9.76 

87 

22.72 

1+2 

9.82 

88 

22.70 

1+3 

9.94 

89 

22.70 

1+1+ 

10.02 

90 

22.69 

1+5 

10.12 
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TABLE  3 

Coefficients  of  the  Two  and  Three 
Parameter  Sections 


Frame 

Two  Parame 

iter  Form 

Three  Parameter  Form 

al 

a3 

al 

a3 

a5 

6l 

0.235739 

0.001735 

0.214840 

0.005425 

0.021767 

91 

0.297847 

-0.075660 

0.292302 

-0.073591 

0.006212 

170 

0.249422 

0.044059 

0.228974 

0.047698 

0.021318 

203 

-0.377228 

0.125725 

-0.462398 

0.059641 

0.107315 

TABLE  4 

Added-Mass  Coefficients  of  four  Sections 
of  a  Cargo  Ship 


Method 

Frame  6l 

Frame  91 

Frame  170 

Frame  203 

C  Cu 

v  H 

C 

V 

CH 

c  c_j 

v  H 

C 

V 

CH 

Two- 

Parameter 

0.997  0.405 

— 

1.139  0.434 

0.937  0.413 

0.777  0.409 

Three- 

Parameter 

0.955  0.487 

1.124 

0 . 466 

0.902  0.489 

0.720 

0.577 

Conformal- 

Mapping 


0.951  0.439 


1.126  0.453 


0.883  0.442 


0 . 7'51  0.611 


IF(DIF.LT.O.O)  GO  TO  16 

15  CONTINUE 

16  DIF  =  ETA(J)  -  W 

AA5  =  B(J)  -  DIF*(B(J)-3(J-1))/(FTA(J)  -  ETA(J-l)) 

T1  =  1.0  +  3.0  *  AA5  *  AA5 
T2  =  3.0  *  TE  +  AA5  *  Tl 

T3  =  2 . 0* ( 1 . 0+CL*TE  )  +  S.O  *  CL  *  SIGMA  /  P 
RAD  =  SQRT(T2*T2  -  4.0  *  Tl  *  T3) 

AL  =  (T2  -  RAD)  /  (2.0  *  Tl) 

AA1  =  (Tl /AL )  -  AA5 
AA3  =  (TE/AL)  -  1.0 
Z  =  3.0  *  (AL  *  AA  5 ) *  *  2 
Z1  =  AL* (AL  -  CL  +  1. )*AA5 
ZZ  =  CL  -  AL  +  1.0 
Z 2  =  2 . 0*  ( ZZ  +  (ZZ-1 . 0)**2 ) 

CV  =  0.5  *  (Z  -  Z1  +  12) 

ONE  =  1.0+  4 . *ZZ  *  ZZ/(3.  *  CL  *  CL) 

TWO  =  ( 5 .  +  40.  *  CL  +  AL*(40. +156. *AA5))*AL*AA 5/(45. *CL*CL) 
CM  =  4.0  *  (ONE  +  TWO)  /  (P*P) 

WRITE(6,17)AA1,AA3,AA5,CV,CH 

17  FORMAT  (  34H  VALUES  OF  A1,A3,A5,CV  AND  CH  ARE/ (  5F  2  0 . 8  )  ) 

CALL  EXIT 

END 
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C  ADDED  MASS  FOR  TWO  DIMENSIONAL  FORMS  BY  CONFORMAL  MAPPING  METHOD 
D I  MENS  ION  A(16),B(lS),CI(16),CF(8),D(8),01{15),F(8),p(g),AA(16), 
1R(8),ANG(91)/R2(91),R0(91),RS(91),XI  (91),ETA(91),Z(91),DZ(91),DRO 
2'(91),TA(91),TE(91),BI(8,9),C(16,16),T(91,16),BB(16,16),CC(16,16), 
3GG (16,16),GA(16,16),P(16),0(16),T1(91,16),AU(91) 

CALL  TRAPS  (-1,-1) 

PP  =  3.14159265 
DANG  =  0.0174533 

RFAD(5,1)((C(I,J),J=1,8),I  =  1,16) 

1  FORMAT ( 8  F 10 . 6  ) 

READ(5,1)((C(!,J),J=9,16)/ 1-1,16) 

READ (5,2)  (<GA(I,J),J  =  1,8), 1=1,16) 

2  FORMAT ( 8  F 10 . 0  ) 

READ ( 5,2  )((GA(I,J),J  =9,16), 1=1,16) 

C  COMPUTATION  OF  COEFF.  B 
READ( 5,149)  AU 

149  FORMAT (40  F2.1) 

100  READ (5, 4)  RS 

4  FORMAT (7F1 1.7) 

BE  =  RS(1) 

HE  =  RS ( 31 ) 

H3  =  DANG  /  3.0 

SO  =  0.0 

DO  77  I  =  1,91 

RS (  I  )  =  RS ( I  )  *  RS (  I  ) 

77  SO  =  SO  +  RS( I )  *  AU ( I ) *  H3 
WRITE (6, 150)  BE  ,  HE  ,  SO 

150  FORMAT (41H  VALUES  OF  BEAM, DRAFT  AND  AREA  OF  SECT  I  ON/ (3F2  0. 8  ) ) 

Z(l)  =  0.0 

DO  7  I  =  1,90 

7  Z( 1+1)  =  Z(l)  +  DANG 
DO  8  I  =  1,91 

8  RS ( I )  =  RS ( I )  *  PP  /  SO 
ES  =  SQRT ( PP/  SO) 

BE  =  BE  *  ES 
HE  =  HE  *  ES 
SO  =  PP 

9  DO  10  J  =  1 >  8 
AJ  =  J 

E  ( J  )  =  0  . 

DO  20  1=1,91 
I  F ( ( J-l ) . ME . 0 )  GO  TO  11 

TE  (  I  )  =-A  LOG  (  SORT  (RS(  I  )  )  )  *  C0S(2.0  *  Z(D) 

GO  TO  20 

11  TE  (  I  )  =((2.*AJ-l.)/(2.*AJ-2.))*OOS(2.*AJ*Z(l  ) ) / RS ( I ) ** ( J-l ) 

20  CONTINUE 

DO  12  L  =  1,90 

E ( J )  =  E ( J )  +  2.*(TE(L)+TE( L+l) )*  DANG 
DO  14  K  =  1,8 
AK  =  K 

TT=  ( 2 . *AJ-1 . ) * ( 2 . *AK-1 . )/ ( 2 . * ( AJ  +  AK-1 .  )) 

B I ( J, K )  =  0.0 


12 
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DO  13  I  =  -1,90 

13  TA ( I )  =  TT*COS ( 2 . * ( AJ-AK) *Z ( I ) )/ RS ( I ) ** ( J+K-l ) 
TA ( 91 )  =  TT*  COS(PP*(AJ-AK))/RS(91)**(J+K-l) 

DO  14  L  =  1,90 

B  !  (  J  ,  K )  =  Bl  (J,K)+2.*(TA(L)+TA(  I.  + 1 )  ) *  D  ANG 

14  CONTI  MU E 
10  CONTINUE 

NEWSH=0 
DO  15  J  =  1,8 

15  B I < J , 9 )  =  E(J) 

DO  18  K  =  2,9 

DA  =  Bl (K-l, K-l ) 

IF (DA. ME. 0.0)  CO  TO  16 

NEWSH-1 

GO  TO  999 

16  DO  17  I  =  K, 9 

17  Bl (K-l, I )  =  Bl (K-l,  I  )  /  DA 
DO  18  J  =  1,8 

I F ( ( K-l ) . EQ. J )  GO  TO  18 
F  =  Bl ( J, K-l ) 

DO  19  I  =  K,  9 

19  Bl (J, I )  =  Bl (J, I )  -  F  *  Bl (K-l, I ) 

18  CONTINUE 

DO  28  1=1,8 
28  B( I )  =  Bl ( I ,9) 

DO  47  I  =  1,91 
R( 1 )  =  SORT ( RS (  I  )) 

CE ( 1 )  =-B ( 1 )  /  R ( 1 ) 

D ( 1 )  =  COS(Z(l  )) 

G ( 1 )  =  S I M ( Z ( I )) 

DO  45  J  =  2,8 

R(  J )  =  R( J  -  1)  *  RS ( I ) 

C F.  ( J  )  =-B  ( J  )  /  R  ( J  ) 

AP  =  2  *  J  -  1 
D  ( J  )  =  COS (AP  *  Z (  I  ) ) 

45  G ( J )  =  S I M ( AP  *  Z ( I ) ) 

X  I  ( I  )  =  R  ( 1 )  *  D  C 1 ) 

ETA( I )  =  R ( 1 )  *  ^(i) 

‘  DO  46  K  =  1,8 

XI ( I )  =  XI ( I )  -  CE ( K )  *  D(K) 

46  ETA ( I )  =  ETA(I)  +  CE(K)  *  G(K) 

RO ( I )  =  SORT (XI(I)*XI(I)  +  ETA ( I )  *  E T A ( I  ) ) 

47  Z (  1  )  =  A TAfi 2  (  ETA(  I  )  , X  I  (  I  )  ) 

SI  =  SO 

DO  48  J  =  1,8 

48  SI  =  SI  +  B ( J )  *  E ( J ) 

ROC12  =  Sl/PP 

RM  =  0.0 

W  =  0.63661977 

DO  49  I  =  1,90 

DZ( I )  =  Z(  1+1)  -  Z (  I  ) 

49.  RM  =  RM+DZO)  *  ( RO  (  I  )  +  RO(l+l))  *  0.5 
RM  =  RM  *  V/ 
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DO  21  I  =  1,8 

21  B  (  I  )  =  B  (  I  )  /  ( RM*  * (  2*1)) 

RM2  =  RM  *  RM 

SI  =  SI  /  RM2 

DO  50  I  =  1,90 

DRO ( I )  =  (ROC  I )  -  RM) /RM 

DO  50  J  =  1,16 

A  ( iJ )  =  0.0 

Q(J)  =  0.0 

P(J)  =  0.0 

H  =  J  *  2 

T( I , J )  =  DRO ( I  )  *  COS (H  *  Z ( I  )  ) 

50  Tl( I ,J)  =  T( I ,J)  *  DRO ( I ) 

DO  51  J  =  1,16 

DO  51  I  =  1,90 

Q(J)  =  Q(  J  )  +  DZ  (  1  )  *  W  *  (TIC  i  ,  <J  )  +  T1  (  I  +1 ,  J  ) ) 

51  P(J)  =  P  (  J  )  +  DZ  (  I  )  *  W  *  ( T  (  I  ,  i! )  +  T(I+1,J)) 

C  COMPUTATION  OF  C(J) 


DO  53  J 

= 

1,16 

SA  =  0.0 

DO  54  K 

s 

1,16 

1 F ( J . EO . 

K) 

GO  TO 

54 

J1  =  ABS ( J 

-  K) 

AL  =  2*K 

- 

1 

SA  =  SA 

+ 

P(K) 

* 

P(J1) 

*  AL 

54 

CONTINUE 

52 

DO  53  K 

= 

1,16 

AJ  =(J  * 

2 

-  1 

)/ 

2 

53 

Cl(d)  * 

-(P(J) 

+ 

SA  - 

Q(J)  * 

S2  =  0.0 

DO  60  J 

= 

1,16 

AJ  =  2  * 

J 

-  1 

60 

S2  =  S2 

+ 

AJ  * 

P(J)  * 

P(J) 

S2  =  SI 

+ 

S2  * 

PP 

*  RM 

*  RM 

ROC22  = 

S2 

/  PP 

DO  62  1 

= 

1,16 

DO  62  J 

= 

1,16 

BB( 1 , J) 

= 

0.0 

62 

GG ( 1 , J  ) 

= 

0.0 

DO  7  3  N 

= 

1,16 

BB(1,M) 

= 

BO!) 

A  ( M  )  =  - 

B  (  N ) 

73 

Dl(l!)  = 

- 

Cl  0! ) 

COM PUT" 

TIOM  OF 

3 

( K  ,  H ) 

DO  63  K  -  2,16 
DO  63  i !  =  1 , 1C 
DO  63  M  =  1,M 
'  L  =  N+l-M 

63  BSC  K , N )  =  SB ( K, N )  +  B ( M )  *  B3(K-1>L) 
C  COMMUTATION  OF  G(J,M) 

DO  64  M  =  2,16 
Ml  =  M  -  1 
DO  64  J  =  1 , Ml 
L  =  M:  -  J 
DO  64  K  =  1 ,  L 
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N  =  L  +  1  -  K 

64  GG ( J , M)  =  GC(J,‘‘)  +  GA(J,K)  *  BS 

C  COMPUTATION  OF  AC  I )  FOR  FIRST 

DO  65  M  =  2,16 
Ml  =  M  -  1 
DO  65  J  =  1,M1 

65  A  ( M )  =  AO-1)  -  A  ( J  )  *  G0( J,M) 

DO  199  I  =  1,16 

DO  199  J  =  1,16 
199  GG( I , J)  =0.0 
C  COMPUTATION  OF  G'(J,M) 

GG ( 1 , 2 )  =  -Cl (1) 

DO  67  M  =  3,16 

N1  =  M  -  1 

FA  =  2  *  Ml  -  1 

G 0 ( N 1 , M )  =  -  EA  *  Cl (1) 

N2  =  M  -  2 
DO  67  J  =  1,M2 
AJ  =  J 

FA  =  2  *  J  -  1 
SF  =  0.0 
L  =  M  -  J 
LI  =  L  -  1 
DO  66  K  =  1 ,  LI 
K1  =  L  -  K 

66  SE  =  SE  +CI (K)  *  Cl C K1 ) 

67  G G ( J , M )  =  FA  * C -C I  CL)  +  AJ  *  SE) 

C  COMPUTATION  OF  D(!) 

DO  GG  M  =  2,16 
Ml  =  r1  -  1 
DO  68  J  =  1 , M 1 

68  Dl(M)  =  D 1 ( M )  -  D1(J)  *  GG(J,M) 

C  COMPUTATION  OF  CORRECTED  A 

DO  59  M  =  1,16 

69  AA(M)  =  A  O')  +  D1(M) 

DO  70  M  =  2,16 

Ml  -  M-l 
DO  70  J  =  1 , M 1 

70  A A ( M )  =  AA(M)  -  (AA(J)  -  D1  ( J  )  ) 
AH  =  0.0 

CV  =  ( 2 . C  *  PM2* ( ROC2 2  +  AA( 1 ) ) 
AA  ( 1 )  =  AA(1)  -  ROC 2 2 
DO  71  I  =  1,16 
DO  71  J  =  1,16 

71  AH  =  AM+  C ( I , J ) * AA ( I )  *  A\(J)  / 
CM  =  16.0  *  PM 2  *  AH  /  (PP  *  HE) 
WRITE (6, 227)  CH,CV 

227  FORMAT  ( 4H  C'1  =  F20.8/4H  CV=  F20.8 
999  WRITE (6, 228)  UFWSH, K 

228  FORMAT (1H  ,215) 

CALL  EXIT 

END 


( K,  M ) 
MAPPING 


*  GG ( J , M ) 

-  1.0)  /  3C**2 

ROC22  **(I+J-1) 
**  2 

) 


FIG.  !a,b  Shapes  of  actual  and  approximate  ship  sections 


71.28 


FIG  lc,d  Shapes  of  actual  and  approximate  ship  sections 
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